Deformation of an asymmetric thin film 
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Experiments have investigated shape changes of polymer films induced by asymmetric sweUing by 
a chemical vapor. Inspired by recent work on the shaping of elastic sheets by non-Euclidean metrics 
[Y. Klein, E. Efrati, and E. Sharon, Science 315, 1116 (2007)], we represent the effect of chemical 
vapors by a change in the target metric tensor. In this problem, unlike that earlier work, the target 
metric is asymmetric between the two sides of the film. Changing this metric induces a curvature 
of the film, which may be curvature into a partial cylinder or a partial sphere. We calculate the 
elastic energy for each of these shapes, and show that the sphere is favored for films smaller than a 
critical size, which depends on the film thickness, while the cylinder is favored for larger films. 

PACS numbers: 46.70.De, 68.60.Bs, 46.25.Cc 



I. INTRODUCTION 

Thin films are three dimensional (3D) objects with one 
dimension much smaller than the other two. Such films 
are not always fiat; they can easily form buckled or wrin- 
kled 3D shapes. Over several years, there has been ex- 
tensive theoretical and experimental research to explore 
the mechanisms for shape selection. This research is im- 
portant to understand the formation of biological struc- 
tures, such as leaves and fiowers, in which thin films as- 
sume well-defined shapes with biological functions. It is 
also important for the design of new synthetic materi- 
als, which should spontaneously form desired morpholo- 
gies [IHl. 

Sharon and collaborators have recently developed an 
important theoretical approach for addressing shape se- 
lection in thin elastic sheets [5l-[l0|. In this approach, a 
film is characterized by a "target metric tensor," which 
describes the ideal spacings between points in the film 
that minimize the local energy. Depending on the math- 
ematical properties of this tensor, there may or may not 
be any global shape of the film embedded in 3D Euclidean 
space that achieves the ideal spacings everywhere. If this 
state is not achievable, the film is geometrically frus- 
trated. Its lowest-energy state will then have residual 
local stresses and strains, and will generally be curved in 
a complex way pdl [l^ . Sharon et al. have demonstrated 
this approach experimentally by using thin films of gels, 
which can be expanded locally by adding a nonuniform 
concentration of a dopant. The gels then relax to the 3D 
shape that has been programmed by the dopant concen- 
tration profile, in agreement with the geometric calcula- 
tions. 

One limitation of Sharon's theory is that it assumes 
the films are uniform across their thickness. This limita- 
tion is significant, because many types of films have some 
variation in their elastic properties across their thickness. 
Indeed, this type of variation might provide an addi- 
tional way to design films to form desired structures. The 



purpose of this paper is to generalize the theory to de- 
scribe asymmetric films, with arbitrary variation across 
the thickness. Through this generalization, we show that 
the asymmetry leads to new types of terms in the elas- 
tic free energy of the film, and we investigate how these 
terms change the curvature. 

As a specific example to motivate this study, we con- 
sider the asymmetric swelling of a thin film by absorption 
of a gas or liquid. A schematic view of this problem is 
shown in Fig. [TJ Here, a cross section through the film 
is illustrated by mass points connected by springs, with 
bold lines indicating that a spring is compressed relative 
to its stress-free length. Before swelling, in Fig. [IJa), all 
the springs are at the stress-free length, and the film is 
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FIG. 1. Schematic illustration of the deformation of a thin 
shell due to swelling: (a) before swelling, (b) after swelling 
without any deformation, and (c) deformed shell. 



2 



flat. After swelling, the intrinsic stress-free lengths of the 
springs vary gradually along the thickness direction, as 
shown in Fig. [ijb). Hence, the film will deform to mini- 
mize the energy, as shown in Fig.[IJc). Because the film is 
3D, the actual geometry is more complex than the cross 
section shown in the figure. It is not obvious whether 
the film should deform into a partial cylinder (with mean 
curvature but no Gaussian curvature) or a partial sphere 
(with both mean and Gaussian curvature). In either case, 
some of the springs will be unable to achieve their intrin- 
sic stress-free length. Hence, this is a simple example of 
a geometrically frustrated structure. 

To address this problem, we extend the approach of 
Sharon et al. to consider a target metric that depends 
on position across the thickness of the film. We calcu- 
late the energy of the deformed film in cylindrical and 
spherical geometries, and find there is a critical lateral 
size, which depends on the film thickness and the gra- 
dient of the intrinsic metric. If the lateral size of the 
film is smaller than this critical size, a partial sphere is 
preferred; otherwise, a partial cylinder has lower energy. 

This paper is organized as follows. In Sec.|TTl we briefly 
review the theory of non-Euclidean plates, and use it to 
calculate the energy in the general case where the target 
metric is proportional to the distance from the midplane. 
We then use this general formula to calculate the energy 
for a cylinder in Sec. IIIII and for a sphere in Sec. IIVI 
Finally, in Sec. El we discuss the results of this study and 
compare with previous work on related elastic problems. 



II. THEORY 

In the first part of this section, we briefly review the 
non-Euclidean theory developed by Sharon et al. for the 
special case in which the target metric is uniform across 
the thickness of the film 0, S]- We then introduce our 
calculation for the asymmetric film with a target metric 
tensor that depends on the distance from the midplane 
of the film. 



A. Brief review of non-Euclidean theory 

Following Sharon et al, we define the elastic free en- 
ergy for any configuration of a material through the fol- 
lowing construction. First, we construct a coordinate 
system x*, for i — 1, 2, 3, in the material frame. The cur- 
rent 3D position of a point (a;^,a;^,x^) is then given by 
). As a result, the current spacing between 
two nearby material points {x^ , x'^ , x^) and {x^+dx^, x'^ + 
dx^,x^ + dx^) is given by (ds)^ — g.ijdx^dx^ (using 
the Einstein summation convention, where Latin indices 
range from 1 to 3, and Greek indices range from 1 to 2). 
Here, the metric tensor is gij — d^iH- d^jll — diR ■ djR. 
By comparison, if we were to cut out a very small piece of 
material around the point {x^,x^,x'^), and let this small 
piece relax to its lowest-energy state, then the spacing be- 



tween those points would be (c?s)^ = (jijdx^dx^ , where gij 
is the target metric tensor that characterizes the intrinsic, 
ideal spacing between the points. The Green-Lagrangian 
strain tensor is then defined as the difference between the 
current metric and the target metric 

= \{gi3 -gij) ■ (1) 

If the material is isotropic, the most general expression 
for the 3D elastic energy density can be written as 

h = Ia^'^^'s^^suu (2) 

where 

Aijki ^ \gijfi + ^i{f^g^^ + f^g^^), (3) 

are the contravariant components of the three- 
dimensional elasticity tensor in curvilinear coordi- 
nates [l^ . If the target metric gij is independent of the 
coordinate x^ across the thickness of the cell, then the 3D 
elastic energy density can be integrated across thickness 
to construct the 2D elastic energy density of the plate, 

/2 = 1 hdx\ (4) 

This 2D elastic energy density includes a stretching en- 
ergy (proportional to thickness h) and a bending energy 
(proportional to h^). The total elastic energy then be- 
comes 

j f2VW\dx'dx^ ^ j h^\dx^dx'dx\ (5) 

integrated over the whole body. 

Note that the target metric tensor gij represents an 
ideal local configuration of the plate with zero stress and 
zero elastic energy density. There might or might not be 
any global configuration in 3D Euclidean space that has 
this target metric tensor everywhere. // such a configu- 
ration exists, then it is a stress-free reference configura- 
tion, which can be used to formulate Truesdell's hyper- 
elasticity principle [l4| . On the other hand, if such a 
configuration does not exist, then there is no such stress- 
free reference configuration. It is still possible to mini- 
mize the elastic energy to find the equilibrium state of 
the plate, but this equilibrium state must be a frustrated 
state with residual local stress. In this case, the object 
can be called a non-Euclidean plane, whose midplane has 
no immersion with zero stretching in 3D Euclidean space. 

B. Asymmetric film 

In our problem of asymmetric swelling of a thin film, 
as shown in Fig. [TJ the target metric must depend on 
position across the thickness of the film. If the gas vapor 
induces expansion of the polymer film, the expansions 
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at the top and the bottom of the film are different due 
to the different concentrations of the gas. Assuming the 
gas concentration varies linearly across the thickness, the 
target tangent vectors can be written as 

= (p + qx^) , (6) 

and hence the target metric is 

gij = S^ ■ gj = kj {f + 2^x3 + {x^f^ . (7) 

In our notation, all symbols with rings represent quan- 
tities before swelling] in particular, gi are the tangent 
vectors and gij is the metric before swelling. Likewise, 
all symbols with bars represent target quantities after 
swelling. The coefficients p and q define how the gas con- 
centration is distributed; p represents a uniform expan- 
sion or compression, and q represents a linear gradient in 
the expansion factor across the thickness of the film. 

From Eq. ([T]) , the strain tensor is defined as the differ- 
ence between the actual metric gij and the target metric 
gij after swelling. To find the actual metric, we must 
consider a specific configuration of the plate. Following 
the second Kirchhoff-Love assumption, we assume that 
points located on any normal to the midplane in the ini- 
tial state remain on that normal in the deformed state, 
but the distance to the midplane may change. Thus, a 
point {x^jx"^,!^ = 0) on the midplane before swelling 
becomes 'R.{x^,x^,x^ = 0) = Rmid(a;^, 2:^) after swelling, 
and a point {x^,x'^,x'^) off the midplane becomes 

R(a;\x2,x3) =R,„id(a;\a;2)+e(a^^)N(x\x2). (8) 

Here, N is the unit normal vector to the midplane, and 
£,(x^) is the new distance to the midplane along the nor- 
mal. Because of swelling, S,{x'^) is no longer equal to x^ . 
To lowest order for a thin film, ^(a;'^) can be written as 
power series 

1 2 

S,{x^) = mx^ + -m' {x^) . (9) 

If we assume that the local volume everywhere remains 
constant during deformation, equal to the swelled local 
volume of the target metric, then the coefficients in this 
expansion are constrained to be m = p and m' = 3q. 

Using this configuration of the plate, we can calcu- 
late the metric tensor gij as a power series in x^. From 
Eqs. ([5]) and the first 2x2 components of gij can be 
written as 

gai3 = aap -2mbai3X^ - m'bap (x^y , (10) 

where a^^ = daJi ■ dpR is the first fundamental form 
of the midplane (i. e. the 2D metric tensor) , and bap — 
N ■ dadpH is the second fundamental form of the mid- 
plane (i. e. the 2D curvature tensor). The first 2x2 com- 
penents of the strain tensor are then Sap = \ {gap — gap)- 
We can now calculate the elastic energy density of the 
plate. Based on the first Kirchhoff-Love assumption that 



the stress is in the local midplane [7|, ll5| , we can express 
the full 3D elastic energy density of Eq. ^ in terms of 
the first 2x2 compenents of the strain tensor as 

/3 = \A''^^'eape^s. (11) 

where 

A''^'^' = 2^i (^-^^^r^-g-^' + r^'f'^ ■ (12) 

If we assume that the local volume everywhere remains 
constant during deformation, then g^-'Sij = 0, and hence 
the elastic modulus A — > cx). As a result, the elasticity 
tensor becomes 

_4a/37<5 ^ 2pL {g°''^f^ + g'^^'f^) . (13) 

Note that the contravariant components of the target 
metric tensor in this expression are the inverse of the 
covariant components of Eq. ([7]). Hence, the elasticity 
tensor depends on x^ as 

pi \^ p pZ \ ' J 

Thus, the 2D elastic energy density can be calculated by 
integrating the 3D elastic energy density over the film 
thickness w, and neglecting terms that has higher order 
than q^ , to obtain 

/w/2 
-w/2 

qw'^ 2D 7 

+ -j^^apOlS + -^aapb^s 

, <1^W^ 2D , l^w^ \ 
H ^2~^a/3^7'5 + 3p2 ^ajSa^S I ■ (i-O) 

To interpret this 2D elastic energy density /2, note 
that it depends on three tensors characterizing the local 
geometry of the midplane: the 2D strain tensor 

= \ - ^ap) , (16) 

which gives the difference between the actual metric and 
the target metric on the midplane, as well as the 2D met- 
ric tensor Oap and the 2D curvature tensor bap. If there 
is no swelling, so that p = I and q = 0, then /2 becomes 
the deformation energy of a symmetric film. In that case, 
the first and second terms in Eq. (ITSt are the stretching 
and bending terms, and the other terms vanish. How- 
ever, if the film is swollen asymmetrically, with g 7^ 0, 
then the last four terms provide new couplings that are 
permitted by the asymmetry. Two of these terms are 
odd in the curvature tensor bap, so they favor sponta- 
neous curvature of the film, which is then coupled to the 
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and 



FIG. 2. (Color online) Illustration of the curved film shapes 
considered in Sections IIIII and IIVI (a) Partial cylinder, 
(b) Partial sphere. 



in-plane strain e^^. The last two terms might seem to be 
higher order in q than the others, but we will show later 
that the equilibrium curvature bap is proportional to q, 
and hence the last five terms are all of order q^. 

For further insight into the 2D elastic energy of 
Eq. ([15]), in the following two sections we use it to calcu- 
late the energy of an asymmetrically swollen film curv- 
ing into a partial cylinder or a partial sphere, as shown 
in Fig. [21 In each case, we calculate the optimal curva- 
ture and determine how it depends on the asymmetry 
parameter q. We then compare the energies of these two 
shapes to see which is favored, as a function of the film 
parameters. 



III. CYLINDER 



We first use Eq. ([T5|) to calculate the energy when the 
film deforms into a partial cylinder, as shown in Fig.[5][a). 
For this calculation, we use an orthogonal curvilinear co- 
ordinate system (^1,^2) on the midplane in the material 
frame, where ^1 is the direction around the circumference 
of the cylinder and ^2 is the direction along the axis. The 
3D position of any point on the midplane can be written 
as 



Rn 



C16 



C26y 



(17) 

where rc is the cylinder radius, Ci and C2 are param- 
eters that measure how much the coordinates elongate 
or shrink, and x, y, and z are unit vectors in Cartesian 
coordinates for the lab frame. From Eq. ([7|), the target 
metric of the midplane is 



Clap 



P 









(18) 



From Eq. I17[ the actual metric and curvature tensors 
(first and second fundamental forms) of the partial cylin- 
der are given by 



CI 






CI 



(19) 



bali — 



( -clh 







(20) 



The elasticity tensor JC^'f^ can be calculated by noticing 
that the target metric tensor before swelling is just the 
identity matrix, independent of position. Hence, the non- 
zero terms in the elasticity tensor are ^^^^^ — ^^2222 _ 
and i"22 ^ _4i2i2 ^ ^42121 ^ ^42211 ^ 2/i. 

We now insert the target metric, actual metric, and 
curvature tensors into the 2D elastic energy of Eq. ([15]) , 
and minimize over the parameters Ci, C2, and r^. To 
lowest order in w, we find that the stretching factors are 
Ci = C2 = p, the cylinder radius is 



3<7' 



(21) 



and the corresponding 2D elastic energy density at that 
minimum is 



29q^_3 



W fJL. 



(22) 



These results imply that the curvature tensor is pro- 
portional to 5, and that the strain tensor = 0; i. e. 
there is no strain on the midplane. Note, however, that 
there is still nonzero shear strain off the midplane be- 
cause the metric tensor off midplane shows anisotropic 
swelling while the target metric favors isotropic swelling; 
that is the reason why the elastic energy is non-zero. In 
the limit of a symmetric film with g ^ 0, then rc — > 00, 
the film remains flat, and there is only uniform swelling 
without energy cost. 

Equation ([^^ shows that the elastic energy density of 
the partial cylinder is uniform in 2D, independent of po- 
sition on the midplane, and hence the total elastic energy 
is just proportional to the film area. In particular, if the 
initial shape of the film is a disk with radius r^ax in the 
material frame, hence radius prmax in the midplane after 
swelling, then the total elastic energy is 



IV. SPHERE 



(23) 



Let us now consider a circular thin film deforming into 
a partial of a sphere, as shown in Fig. [2][b). For this 
problem, it is convenient to use polar coordinates (r, cj)) 
on the midplane in the material frame, where r is the 
radial displacement from the center of the circular film 
before swelling, and is the azimuthal angle, which is 
assumed not to change during swelling and deformation. 
The 3D position of any point on the midplane can then 
be written as 

R-mid — fs sin 9{r) cos ^x-t-r^ sin d{r) sin (j)y + rs cos 0{r)z, 

(24) 
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where is the radius of the sphere and 9{r) is a monoton- 
ically increasing function of r, to be determined, which 
describes how the material stretches or shrinks in the 
radial direction. (In particular, 9(r) gives the angular 
position on the partial sphere up from the (— 2;)-axis cor- 
responding to the radial position r in the material frame.) 
In this coordinate system, the target metric of the middle 
plane is 



the total elastic energy is 





9 9 



(25) 



Note that this target metric is equivalent to Eq. ([T8l) . 
but in a different coordinate system. From Eq. (|24|) . the 
actual metric and curvature tensors (first and second fun- 
damental forms) are given by 



and 



ba/B = 



9'(r)2 












9 • 9 

rz sin 



e{r) 



(26) 





-rs sin^ 9{r) 



(27) 



In this coordinate system, the non-zero components of 
the elasticity tensor ^"'^75 _^iiii ^ ^42222 ^ 
4/x/r4, and ^1^22 ^ _4i2i2 ^ _42i2i ^ _422ii ^ 2^JL/r■'. 

The 2D elastic energy of Eq. ([T5|) is now a functional 
of the stretching function 9{r) and the sphere radius Tj. 
To minimize this energy over 9(r), we expand 9{r) as a 
power series in r and minimize over the series coefficients. 
The leading terms are then 



9{r) 



ip- 



3 2 
p w 

6r? 



8r, 



(28) 



To minimize the energy over r^, we must consider two 
regimes in terms of q, w, and rmax, the radius of the 
film in the material frame. In the first regime, where 
qr^s^^/w <C 1, the optimum sphere radius is — p^/q, 
and the total elastic energy is 



V. DISCUSSION 



(30) 



Our calculations in the preceding sections lead to 
specific conclusions about the cylindrical and spherical 
shapes, as well as more general insight into elastic theory 
for asymmetric films. 

For the specific problem of cylindrical and spherical 
shapes, we can see that Eq. ((23)) for the elastic energy of 
a partial cylinder is between the two regimes of Eqs. ([29]) 
and (pO)) for a partial sphere. This result implies that 
the spherical deformation is favored for disks of small 
radius Tmax < ^critical, whilc the cylindrical deformation 
is favored for disks of large rmax > ''critical- Calculating 
the actual value of rcriticai requires higher-order terms 
than we have presented here, and the numerical result is 
''critical = 1 ■6{wp/ q)^ ^"^ . To Understand this result, note 
that the partial sphere has an isotropic deformation in 
the midplane, which is consistent with the target met- 
ric, while the partial cylinder breaks is anisotropic in the 
midplane and disagrees with the target, which costs ex- 
tra energy. For that reason, the partial sphere is favored 
for small r^ax- By contrast, for large rmax the partial 
sphere must have extra stretching in the midplane, and 
hence it becomes disfavored with respect to the partial 
cylinder. 

More generally, we have shown that the theoretical ap- 
proach of Sharon and collaborators can be applied to 
asymmetric films. Through this approach, we transform 
the 3D elastic energy into the effective 2D elastic en- 
ergy of Eq. p5|) . This effective 2D elastic energy shows 
the standard stretching and bending energies, which have 
been studied extensively for symmetric films, as well as 
new terms arising from the asymmetry. These new terms 
include a spontaneous curvature term, which is linear in 
the curvature tensor and hence favors curvature of the 
asymmetric film, as well as a coupling between sponta- 
neous curvature and in-plane strain. These new terms 
should provide new opportunities to design synthetic ma- 
terials that will spontaneously form desired shapes for 
technological applications. 



(29) 
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